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Abstract. The ground state of a gas of Bosons confined in an external trap 
potential and interacting via repulsive two- body forces has recently been shown 
to exhibit complete Bose-Einstein condensation in the dilute limit, yielding 
for the first time a rigorous proof of this phenomenon in a physically realistic 
setting. We give here an account of this work about the Gross-Pitaevskii limit 
where the particle number N goes to infinity with Na fixed, where a is the 
scattering length of the interparticle interaction, measured in units of the trap 



1. Introduction and Main Results 

During the last few years it has become experimentally feasible to realize the 
long-predicted Bose-Einstein condensation (BEC) of gases by confining them in 
traps at very low temperatures. A rigorous theoretical demonstration of this phe- 
nomenon ~ starting from the basic many-body Hamiltonian of interacting particles 



- is, however, a very difficult task. Following [LS|, we will provide in this paper 
such a rigorous derivation for the ground state of Bosons in an external trap in- 
teracting with repulsive pair potentials, and in the well-defined limit in which the 
Gross-Pitaevskii (GP) formula is applicable. It is the first proof of BEC for inter- 
acting particles in a continuum (as distinct from lattice) model and in a physically 



realistic situation. We present here a detailed version of the proof given in [LS 



We remark that an extension of the results presented here was recently obtained 



in |LSY4], where it was shown that the ground state is 100% superfluid in the GP 
limit. 

The Gross-Pitaevskii limit under discussion here is a mathematically simpler 
limit than the usual thermodynamic limit in which the average density is held fixed 
as the particle number goes to infinity. A proof of BEC in this limit is still missing. 
The only available rigorous results concern the non-interacting gas, various mean- 



lattice gas |KLS] 



field and other toy models (see, e.g., |LVZ] for a recent preprint), or the hard core 
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In the GP limit one also lets the range of the potential go to zero as N goes 
to infinity, but in such a way that the overall effect is non-trivial. That is, the 
combined effect of the infinite particle limit and the zero range limit is such as to 
leave a measurable residue — the GP function. The limit in which the GP function 
can be expected to be equal to the condensate wave function should be chosen so 
that all three terms in the GP functional ( |l.2D make a contribution. This indicates 
that fixing Na as — > oo is the right thing to do, and this is quite relevant since 
experimentally N can be quite large, 10® and more, and Na can range from about 



1 to 10^ |PeSm]. Fixing Na also means that we really are dealing with a dilute 
system, because the mean density p is then of the order N and hence a^p ^ N~'^, 
meaning that the range of the interaction is much shorter than the mean particle 
distance. 

We shall now describe the setting more precisely. Let H be the Hamiltonian for 
identical Bosons in a trap potential V in R^, interacting via a pair potential v: 

N 

i/ = 5] ( - A, + y(x,)) + ^i^^-^,)- (1-1) 

1=1 l<i<j<N 

It acts on the subspace of totally symmetric functions in {^^L^(K'^). Units are 
chosen such that h = 2m = 1, where m denotes the particle mass. We assume the 
trap potential F to be a locally bounded function that tends to infinity as |x| oo. 
The interaction potential v is assumed to be nonnegative, spherically symmetric, 
and have a finite range Rq. Its scattering length will be denoted by a. (For the 
definition of scattering length, see [LSY1| or | LY2| ].) The finiteness of the range 



will be assumed for simplicity, but this assumption can be relaxed. Note that we 
do not demand v to be locally integrable; it is allowed to have a hard core, which 
forces the wave functions to vanish whenever two particles are close together. In 
the following, we want to let a vary with TV, and we do this by scaling, i.e., we 
write t'(x) = z)i(x/a)/a^, where vi has scattering length 1, and keep vi fixed when 
varying a. The ground state energy of H, denoted by E'-^^, then depends only on 
A'' and a (for fixed V and vi), so the notation E'^^{N,a) is justified. The results 
below are independent of the particular shape of the interaction potential vi . 
The Gross-Pitaevskii functional is given by 



^^''M - J (|V0(x)p + y(x)|(/)(x)|2 + .9|0(x)|4)d3x . (1.2) 
Here g is a positive parameter that is related to the particle number and the scat- 



tering length of the interaction potential appearing in (1.1) via 

g^A-KNa. (1.3) 

We denote by (j>'^^ the minimizer of S'^^ under the normalization condition J |</)p = 
1. Existence, uniqueness, and some regularity properties of were proved in the 



appendix of [ LSYl j. In particular, (j)'^^ is continuously differentiable and strictly 
positive. Of course cj)'^^ depends on g, but we omit this dependence for simplicity 
of notation. For later use, we define the projector 

P«P = |O(0^"|. (1.4) 
The minimizer fulfills the variational equation 

-A0«P + y</>«P -I- 25(0«P)^ = A^GP^GP ^ (^ 5) 
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which is called the GP equation. Here /i^^ is the chemical potential, given by 

f,^^ = E^^+g [id^^^W (1.6) 



with = E^^{g) the lowest energy of f under the condition J = 1. 

It was shown in ]LSY1 | (see Theorem ^ below) that, for each fixed g, the 
minimization of the GP functional correctly reproduces the large N asymptotics 
of the ground state energy and density oi H ~ but no assertion about BEG in this 
limit was made. This was extended in |LS|, where complete BEG was proved. The 
precise statement is as follows. 

Let ^' denote the nonnegative and normalized ground state of H. BEG refers 
to the reduced one-particle density matrix 

7(x, j ^^(x, X)vl'(x', X)dX , (1.7) 

where we denoted X = (x2, . . . , x^v) and dX = Y\^^2 '^^^j ^'^^ short. It is the kernel 
of a positive trace class operator 7 on L^(M'^), with Tr[7] = N. 

Gomplete (or 100%) BEG is defined to be the property that ■^7(x, x') becomes 
a simple product (p{x.)(p{x') as iV — > 00, in which case ip is called the condensate 
wave function. I.e., -^-7 converges to a one-dimensional projection. In particular, 
there is one eigenvalue of 7 of the order N (even equal to N in the limit), and all 
the others are of lower order. In the GP limit, i.e., N ^ 00 with g — AnNa fixed, 
we can show that this is the case, and the condensate wave function is, in fact, the 
GP minimizer cj)'^^ . 

Theorem 1 (Bose-Einstein Condensation) . For each fixed g 

lim l7(x,x') = 0«P(x)0GP(x') (1.8) 

in trace class norm, i.e., Tr[|7/iV-FGP|] -> 0. 

Remark 1. Theorem |l| implies that there is 100% condensation for all n- 
particle reduced density matrices of 4", i.e., they converge to the one-dimensional 
projector onto the corresponding n-fold product of (f)'^^ . Sec |LS| for details. 

Remark 2. Gonverg ence does not hold for the Sobolev norm Tr[|(l — A)(7/iV — 
P°P)|]. It can be shown @ that 

hm ^Tr[-A7] = Tr[-ApGP] + .gs / (1.9) 

for some parameter < s < 1 depending only on the interaction potential vi . 

A corollary of Theorem |l|, important for the interpretation of experiments, 
concerns the momentum distribution of the ground state. 

Corollary 1 (Convergence of momentum distribution). Let 

p{k) = J 7(x, x') exp[ik • (x - xOld^xd^x' (1.10) 
denote the one-particle momentum density of . Then, for fixed g, 

hm lp(k) = |0«P(k)p (1.11) 

N^ca IS 

strongly in L^{M.^). Here, (fP^ denotes the Fourier transform of (fP^ . 
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Proof. If T denotes the (unitary) operator 'Fourier transform' and if h is an 
arbitrary L°°-function, then 

^Jph-J I^^^P/^ = \TT[T-\^/N^P''^)Th]\ < ||/.||ooTr[|7/^-pGP|] ^ 

(1.12) 

from which we conclude that 

\\p/N - |0GP|2||^ < Tr[|7/iV - PGP|] . (1.13) 

□ 

It is important to note that in the limit considered v becomes a hard potential 
of short range. This is the opposite of the usual mean field limit, where the strength 
of the potential goes to zero while its range tends to infinity. 

We also wish to emphasize that in this GP limit the fact that there is 100% 
condensation does not mean that no significant interactions occur. The kinetic and 
potential energies can differ markedly from those obtained with a simple variational 
function that is an A^-fold product of one-body condensate wave functions. This 
assertion might seem paradoxical, and the explanation is that near the GP limit 
the region in which the wave function differs from the condensate function has a 
tiny volume that goes to zero as N ~^ oo. Nevertheless, the interaction energy, 
which is proportional to iV, resides in this tiny volume. 

Before proving Theorem |^, let us state some prior results on which we shall 
build. These results concern the asymptotic behavior of the ground state energy 



and density of (IT) in th e limit N ^ oo with g = AirNa fixed. The following 



Theorem ^ was proved in ||LSY1 



Theorem 2 (Asymptotics of Energy and Density). Let p{x) — 7(x, x) 
denote the density of the ground state of H . For fixed g = AnNa, 

hm li?Q^(^,a)=i?«P(5) (1.14) 

N^cxs I\ 

and 

hm lp(x) = |0«P(x)p (1.15) 
in the sense of weak convergence in 



Remark 3. The convergence in ( |1.15 ) was shown in | LSYl[ | to be in the weak 



X"'^(M^) sense, but our result here implies strong convergence, in fact. This can 
easily be deduced from Theorem |l| (cf. the proof of Corollary 

We now give a brief outline of the proof of Theorem |l|. There are two essential 
ingredients. The first is a proof that the part of the kinetic energy that is associated 
with the interaction v is mostly located in small balls surrounding each particle. 
More precisely, these balls can be taken to have radius A^~^/^^, which is much 
smaller than the mean particle spacing N"^^^. This allows us to conclude that the 
function of x defined for each fixed value of X by 

1 

(/)GP(x) 

has the property that Vx/x(x) is almost zero outside the small balls centered at 
points of X. This is made precise in Lemma |l|. 



/x(x) = -gp7::t*(^'X)>o (1.16) 
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The complement ol the small balls has a large volume but it can be a weird set; 
it need not even be connected. Therefore, the smallness of Vx/x(x) in this set does 
not guarantee that /x(x) is nearly constant (in x), or even that it is continuous. 
We need /x(x) to be nearly constant in order to conclude BEC. What saves the day 
is the knowledge that the total kinetic energy of /x(x) (including the balls) is not 
huge. The result that allows us to combine these two pieces of information in order 
to deduce the almost constancy of /x(x) is the generalized Poincare inequality in 
Lemma 0. 



2. The Proof 



We start with the following Lemma. It shows that to leading order all the 
interaction energy is concentrated in small balls. The proof we give here is a 
detailed version of the one sketched in [LS|. Instead of referring to the methods of 



LSY1|, we follow a slightly different route and use the (partly simpler) methods 



of the review article [ LSSoY | . 

Throughout the paper we suppress the dependence on N for simplicity of no- 
tation. For instance, in Eq. (2^) both X and /x depend on N. 



Lemma 1 (Localization of the Energy). For fixed X. let 



rJx = <^ X G 



Then 



minlx-Xfcl > iV"^/"l 

k>2 J 



hm fdxf d3x|0GP(x)p|Vx/x(x)P =0 



(2.1) 



(2.2) 



Note that ( p.2[ ) docs not imply that the kinetic energy of /x goes to zero. In 
fact, it can be shown |LS, CS| that if one replaces fix by M"^ in (2.2) the result is 



hm fdxf d^x|0«P(x)|2|Vx/x(x)| 



gs I |</.«P(x)|4d3x 



(2.3) 



(compare with (1.9)). Here < s < 1 is a parameter depending only on the 
interaction potential v. (For example, s = 1 in the case of hard core Bosons; in 
general, s < 1.) The right side of ( ^.3D is 0(1) in the limit considered, and is not 
necessarily small. Thus Lemma ^ shows that all the interaction energy is localized 
in small balls of radius N^'^^^'^ surrounding each particle. 



Proof of Lemma |T[ We shall show that 
fdxf d3x|</.GP(x)nVx/x(x)P 



X - Xfcl 



25l0^^(x)| 



fc>2 

> -9 



4.GP/ 



(x)l^d-^x-o(l) (2.4) 



as ^ oo. Here il^ the complement of fix in R"^. We claim that this implies 
the assertion of the Lemma. To see this, note that the left side of (2.4) can be 
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written as 



±EQM _ GP 

N ^ 



[dxf d3x|0GP(x)|2|v^/x(x)|2 , (2.5) 



where we used partial integration and the GP equation (L5), and als o th e symmetry 
of The convergence of the energies in Theorem ^ and the relation (1^) now imply 

the desired result. 

We are left with the proof of (2^). This is actuall y just a deta iled examination 
of the lower bounds to the energy derived in [LSY1| and [ LY1|, and we use the 
methods from there. The proof presented here follows closely [LSSoY and differs 



slightly from |LSY1 



Writing /x(x) = Ilk>24''~'^ {^k)F{x, X) and using that F is symmetric in the 
particle coordinates, we see that (^^) is equivalent to 



^Q(^) > -9 



/GP|4 



r-o(i) 



(2.6) 



where Q is the quadratic form 

TV N 

QiF)^J2 / jV.Fp[]|</.GP(xfe)pd^x, 



i=l 
N 

E 

1=1 



^t;(x,-x,)-25|0«P(xO| = 



N 



\F\'Yl\q^^^i^,)\^d'^,, (2.7) 



k=l 



with = {{xi^X) G M^wi jninfc-,, |x, - xfcl < N-y^^}. 

While (2.6) is not true for all conceivable -F's satisfying the normalization con- 
dition 



/ \F\^l[\<^^^i^,)\'d'^,^l, 
fe=i 



(2.8) 



it is true for an F, such as ours, that has bounded kinetic energy (2.3). In fact, we 
will show that 



1 

TV 



Af . JV 

Q{F)+ey2 |V,i^pTT|0GP(x,)|2d3^, 

k=l 



>-9 -0(1) , (2.9) 



for all F satisfying (2.8), with e = o(l) as — > oo. 

To estimate the left side of (2.9) from below, we divide space into boxes of 
side length L, labeled by a, and distribute the particles over the boxes. If we use 
Neumann boundary conditions in each box, this can only lower the energy, since 
we effectively allow discontinuous functions. Moreover, since v is positive, we can 
neglect interactions am ong particles in different boxes for the lower bound. To be 
precise, the left side of (2^) is bounded below by 



1 inf yinf^^" 



(2.10) 
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where the infimum is over all distributions of Ua particles in the boxes a, under the 
constraint that '^a — -^j ^^'^ 



E 

i=l 



|V,F,| 



|V,F„| 



fc=i 



+ 



aGP/ 



(xfe)rd"^Xfc . (2.11) 



k=l 



Here all the integrals are restricted to the box a, and Fa is a function ofua variables. 

We now fix some M > 0, that will eventually tend to oo, and restrict our- 
selves to boxes a inside a cube Am of side length M. Since v > the con- 
tribution to (2.10) of boxes outside this cube is easily estimated from below by 
— 2(7A^suPx^Am l</'*^^(x)P, which, divided by N, is arbitrarily small for M large, 
since ^'^^ decreases faster than exponentially at infinity (|LSY1], Lemma A. 5). 

For the boxes insi de the cube Am we want to use the results on the homogeneous 
Bose gas obtained in [ LYl | , and therefore we must approximate by constants 
in each box. Let /Omax and Pmin, respectively, denote the maximal and minimal 
values of \(f>' 



GP|2 



in box a. Define 

^'a(xi, . . . ,X„^) = Fa{xi, 



,x„jn<^«p(x,) 



k=l 



and 



*«(X1, . . . ,X„J = ^^„(X1, . . . ,X„J [] </.«P(Xfc) 



(2.12) 



(2.13) 



fe=i 



We have, for all 1 < i < 



V,.VI,W|2 



ni0«p(x,)pd3^. 



anw; 



IT. 



k=l 



v{x. 



J] d^Xk . (2.14) 



fc=i 



We now use the following Lemma, which was proved in [ LYl | . It allows us to replace 
w by a 'soft' potential, at the cost of sacrificing kinetic energy and increasing the 
effective range. 

Lemma 2. Let v{x) > with finite range Ro- Let U{r) > be any function 
satisfying JU{r)r^dr < 1 and Uir) — for r < Rq. Let B C M.^ be star shaped 
with respect to (e.g. convex with E B). Then for all functions ip 



(2.15) 



By dividing a for given points Xi , . . . , x„^ into Voronoi cells Bi that contain 
all points closer to x^ than to Xj with j ^ i (these cells are star shaped w.r.t. x^, 
indeed convex), and choosing a U with radius R < TV^''/^'', we see that (2.14) is 
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bounded below by 



J a Pmax J a 



(2.16) 



where ti is the distance of x.j to its nearest neighbor among the other points Xj, 



ii(xi, . 



, ,x„^ j = mm |Xi - Xj 



As in [LYl] we choose for U the potential 

'3(i?3 - i?3)-i fori?o<r<i? 
otherwise , 



Uir) = 



(2.17) 



(2.18) 



with R determined by Qx SiS R = N Note that R^ Rq for N large enough, 

since i?o = 0(-/V^^) in the limit considered. 
Since ^'q — 0'^^(xi)^'a'* we can estimate 

iV.^-cl' < 2p„,ax|V,*«|2 + 2|*«|2Cm (2.19) 

with 

Cm = sup |V(/.°P(x)|2 . (2.20) 

xGAm 



Inserting (2.19) into (2.11), summing over i and using p^^(xi) < pmax in the last 
term of ( 2.11 ) , we get 



Qa{Pa) ^ Pmin j^(J 



\F„ 



where {ua, L) is the ground state energy of 

£ ( - ieA, + aUiU)) 

i=l 



(2.21) 



(2.22) 



in a box of side length L. We want to minimize (2.21) with respect to tIq,, and drop 
the subsidiary condition X]q = N in (2.10). This can only lower the minimum. 
For the time being we also ignore the last term in ( 2.21 ). The total contribution 
of this term for all boxes is bounded by eCmN and will turn out to be negligible 
compared to the other terms. 



It was shown in |LY1] (see also |LSSoY]) that 



(2.23) 



with Yq, = a^najl? 1 provided e > YJ and Ua > (const. )¥« • The condition 
on e is certainly fulfilled if we choose e — Y^l^"^ with Y ~ a^N/L^. We now want 
to show that the minimizing the right side of (2.21) is large enough for (2.23) 
to apply. 

If the minimum of the right side of (2.21) (without the last term) is taken for 
some fia^ we have 

^ {EY{n„ + 1, i) - ("a, i)) > 25P.nax ■ (2.24) 

Pmax 

On the other hand, we claim that 
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Lemma 3. For any n 



E^{n + 1, L) - E^{ti, L) < Sna^ . (2.25) 



Proof. Denote the operator (2.22) by i/„, with Ua = n, and let be its 
ground state. Let t'^ be the distance to the nearest neighbor of among the n + 1 
points Xi, . . . , Xn+i (without x^) and ti the corresponding distance excluding x„+i. 
Obviously, for 1 < i < n, 

U{t',) <UiU) + Ui\x,-Xn+i\) (2.26) 

and 

n 

f/(C+i)<E^(l^'-^"+il) ■ (2.27) 
1=1 

Therefore 

n 

Hn+1 <Hn- ieA„+i + 2aJ2 U{\^^ - x„+i |) . (2.28) 

i=l 

Using '^n/L^/'^ as trial function for Hn+i we arrive at ( ^.21^ ). □ 

Eq. ( 2.25| ) together with ( ^.24| ) shows that is at least '-^ NpniaxL^ (recall 
that g = AnNa). We shall choose L ^ iV^^/^", so the conditions needed for ( ^.23 ) 
are fulfilled for N large enough, since Pmax = 0(1) and hence fia ^ N'^/^^ and 

In order to obtain a lower bound on (2.21), we can use Ya < Y in the error 
term in (2.23). We therefore have to minimize 

4na (i _ CY^nA - 2n^Np„,,^ , (2.29) 

and we can drop the requirement that has to be an integer. The minimum of 
( p.29| ) is obtained for 

Using ( |2l0| ), ( p^ and ( |2^ ) this gives the following lower bound on the left side 
of (2.£), including now the last term in (2.21) as well as the contributions from the 
boxes outside A^/: 

-9 E ^-i„i' f4^7r^^7T7m ) - 2.9 ^up |</.°P(x)|2 - eCu ■ (2.31) 

acAif \ftninU-'--J^ J/ x^A^, 

Now is differentiable and strictly positive. Since all the boxes are in the fixed 
cube Am there are constants C" < oo, C" > 0, such that 

Pmax ^ Pmin ^ C L, Pmin ^ C . (2.32) 

Since L ^ iV-Vio ^nd Y - A^-i'^/io we therefore have, for large N, 

(2.33) 



Also, 

9 E pLnL^<gl\c^^^\'. (2.34) 



aCAjv 
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Hence, for some constant depending only on g and M , 

(PI > -g I |(/)^P|*-(const.)iV"'/'"-2(7 sup |0°P(x)p 
This proves the desired result, and finishes the proof of Lemma |^. 



(2.35) 

□ 



In the following, IC C denotes a bounded and connected set that is suffi- 
ciently nice so that the Poincare-Sobolev inequality (see [ LLo , Thm . 8.12 ]) holds 
on /C. In particular, this is the case if K, satisfies the cone property [LLo| (e.g., if 
/C is a ball or a cube). The following Lemma is a generalization of the Poincare 
inequality. It can be further generalized to the case, and, with a different and 
more complicated proof, to the case of magnetic fields |LSY3|. 

Lemma 4 (Generalized Poincare Inequality). For m > 2 let fZ C K™ he 

as explained above, and let h be a bounded function with J/^h — 1 . There exists 
a constant C ( depending only on K, and h ) such that for all sets G K. and all 
f G H^{IC) with Jj^ fhd™:x. = 0, the inequality 



JjVf{K)fd^K+ (^l^y^'"|^|v/(x)|2d'"x> i^|/(x)|2d'"x (2, 



36) 



holds. Here \ ■ \ is the volume of a set, and D,'^ = JC\Cl. 



Proof. By the usual Poincare-Sobolev inequality on IC (see [LLo, Thm. 8.12]), 



WfWhiK) < C||V/||i.™/,„^.)(^) (2.37) 

for some constant C, if m > 2 and //i = 0. Using the triangle inequality we can 
estimate 

WIWUk) < 2C (||V/||i.„/(„+.,(^) + I|V/||i.„/<„+„(f,.)) . (2.38) 
Applying Holder's inequality 

||V/L.w(™+.)(o) < ||V/|U2(o)|r!r/™ (2.39) 

(and the analogue with fl replaced by fl'^), we see that ( |2.36D holds with C 
2|/C|2/'"C. □ 

The important point in Lemma ^ is that there is no restriction on ft concerning 
regularity or connectivity. Combining the results of Lemmas |l] and ^ we now are 
able to prove Theorem |l|. 

Proof of Theorem |i|. For some M > let AC = {x e R^, [x] < M}, and 
define 

(/x). - ;^|^cP;.)p,3. I |^-(x)|Vx(x)d3x . (2.40) 

We shall use Lemma |, with m = 3, /i(x) = |0C!P(x)|2/ |0C!P|2^ fl = Qy^ n K. and 
/(x) = /x(x) — {fx.)ic (see (2.1) and (1.16)). Since (p'^^ is bounded on IC above and 
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below by some positive constants, this Lemma also holds (with a different constant 
C") with d^x replaced by \(f)'^^ {x.)\^d^x in (0|). Therefore, 



J dX^d3^|0GP(x)p[/x(x) - (/x);c]' 



<C" dX 



°P(x)P|V./x(x)pd3x 



M2 



'(x)nVx/x(x)|^d3x 



, (2.41) 



where we used that n /C| < (47r/3)A^^^/^^. The first integral on the right side 
of ( 2.4l| ) tends to zero as — > oo by Lemma |^, and the second is bounded by (2.3). 
We conclude that 



hm fdxf d3x|0GP(x)|2[/x(x) - (/x);c]' - 

Moreover, since 

/ |0GP(x)p/x(x)d3x< / |0GP(x)p/x(x)d3x 
by the positivity of /x, 



/ |0«P(x)pd3x 



Hence, by ( ^l2| ). 



1 



liminf-(0GP| uGP 



)> / |0°^(x)|2d3x lim / dX / d3x|^'(x,X)| 



(2.42) 



(2.43) 



(2.44) 



(2.45) 



It follows from ( 1.15| ) that the right side of this inequality equals [/^ (x) pd'^x] 
Since the radius of JC was arbitrary, we conclude that 



lim A(^GP|^|^GP^^i 

A'— s-oo I\ 



(2.46) 



implying convergence of '^/N to P'^^ in Hilbert-Schmidt norm. Since the traces 
are equal, convergence even holds in trace class norm (cf. |S^, Thm. 2.20]), and 
Theorem ^ is proven. □ 

Throughout the paper we were dealing with Bosons in three-dimensional space. 
However, the method presented here also works in the case of a 2D Bose gas. 
The relevant parameter to be kept fixed in the GP limit is 5 = ATrN/\\n{a'^N)\, 
all other considerations carry over without essential change, using the results in 



LSY2, LY2|. We also point out that our method necessarily fails for the ID Bose 



gas, where there is no BEG in the ground state [PiSt|. An analogue of Lemma 1 
cannot hold in the ID case since even a hard core potential with arbitrarily small 
range produces an interaction energy that is not localized on scales smaller than 
the mean particle spacing. 
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